Theory of acoustic surface plasmons 
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Recently, a novel low-energy collective excitation has been predicted to exist at metal sur- 
faces where a quasi two-dimensional (2D) surface-state band coexists with the underlying three- 
dimensional (3D) continuum. Here we present a model in which the screening of a semiinfinite 3D 
metal is incorporated into the description of electronic excitations in a 2D electron gas through the 
introduction of an effective 2D dielectric function. Our self-consistent calculations of the dynam- 
ical response of the 3D substrate indicate that an acoustic surface plasmon exists for all possible 
locations of the 2D sheet relative to the metal surface. This low-energy excitation, which exhibits 
linear dispersion at low wave vectors, is dictated by the nonlocality of the 3D dynamical response 
providing incomplete screening of the 2D electron-density oscillations. 



PACS numbers: 71.45.Gm, 73.20.At, 73.50.Gr, 78.47.-l-p 
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I. INTRODUCTION 

Since the early suggestion of Pinegi- that low-energy 
plasmons with sound-like long-wavelength dispersion 
could be realized in the collective motion of a system 
of two types of electronic carriers, these modes have 
spurred over the years a remarkable interest and research 
activityi^ 

The possibility of having a longitudinal acoustic mode 
in a metal-insulator-semiconductor (MIS) structure was 
anticipated by Chapliki^ Chaplik considered a simpli- 
fied model in which a two-dimensional (2D) electron gas 
is separated from a semiinfinite three-dimensional (3D) 
metal. He found that the screening of valence electrons in 
the metal changes the 2D plasmon energy from its char- 
acteristic square-root wave- vector dependence to a linear 
dispersion, which was also discussed by Gumhalter— in his 
study of transient interactions of surface-state electron- 
hole (e-h) pairs at surfaces. 

Nevertheless, acoustic plasmons were only expected to 
exist for spatially separated plasmas, as pointed out by 
Das Sarma and Madhukari^ The experimental realiza- 
tion of two-dimensionally confined and spatially sepa- 
rated multicomponent structures, such as quantum wells 
and heterostructures, provided suitable solid-state sys- 
tems for the observation of acoustic plasmonsi^ Acoustic 
plasma oscillations were then proposed as possible candi- 
dates to mediate the attractive interaction leading to the 
formation of Cooper pairs in high-Tc superconductors.'''^ 

Recently, Silkin et al^ have shown that metal surfaces 
where a partially occupied quasi-2D surface-state band 
coexists in the same region of space with the underly- 
ing 3D continuum support a well-defined acoustic surface 
plasmon, which could not be explained within the orig- 



inal local model of Chapliki^ This low-energy collective 
excitation exhibits linear dispersion at low wave vectors, 
and might therefore affect e-h and phonon dynamics near 
the Fermi level »iS 

In this paper, we present a model in which the screen- 
ing of a semiinfinite 3D metal is incorporated into the 
description of electronic excitations in a 2D electron gas 
through the introduction of an effective 2D dielectric 
function. We find that the dynamical screening of va- 
lence electrons in the metal changes the 2D plasmon en- 
ergy from its characteristic square-root behaviour to a 
linear dispersion, not only in the case of a 2D sheet spa- 
tially separated from the semiinfinite metal, as antici- 
pated by Chaplik,^ but also when the 2D sheet coexists 
in the same region of space with the underlying metal, as 
occurs in the real situation of surface states at a metal 
surface. Furthermore, our results indicate that it is the 
nonlocality of the 3D dynamical response which allows 
the formation of 2D electron-density acoustic oscillations 
at metal surfaces, since these oscillations would otherwise 
be completely screened by the surrounding 3D substrate. 

Unless stated otherwise, atomic units are used 
throughout, i.e., = h = me = 1. 



II. THEORY 

A variety of metal surfaces, such as Be(OOOl) and the 
(111) surfaces of the noble metals Cu, Ag, and Au, are 
known to support a partially occupied band of Shockley 
surface states with energies near the Fermi level. ■'^'^ Since 
the wavefunction of these states is strongly localized near 
the surface and decays exponentially into the solid, they 
can be considered to form a 2D electron gas. 
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In order to describe the electronic excitations occurring 
within a surface-state band that is coupled with the un- 
derlying continuum of valence electrons in the metal, we 
consider a model in which surface-state electrons com- 
prise a 2D electron gas at z = Zd {z denotes the coor- 
dinate normal to the surface), while all other states of 
the metal comprise a 3D substrate consisting of a fixed 
uniform positive background (jellium) of density 



n+{z) 



n, z < 



0, elsewere. 



(1) 



plus a neutralizing inhomogeneous cloud of interacting 
electrons. The positive-background charge density n is 
often expressed in terms of the 3D Wigner radius r^^ — 
(3/47rn)^/3/ao, ao = 0.529A being the Bohr radius. 

We consider the response of the interacting 2D and 3D 
electronic subsystems to an external potential (^"^^'(r, ui). 
According to time-dependent perturbation theory, keep- 
ing only terms of first order in the external perturbation, 
and Fourier transforming in two directions, the electron 
densities induced in the 2D and 3D subsystems are found 
to be 



(jy^^'-^^z-^q^bj) + I dz'v{z,z';q)Sn3D{z';q,i^) 



and 



Sn3D{z;q,uj) = J dz' X3D{z,z';q,uj) 
<j,--\z'- q,u;)+ [ dz"v{z', z"; q)5n2D{z"; q, uj) 



(3) 



Here, q is the magnitude of the 2D wave vector parallel 
to the surface, X2d(9,w) and xsoiz, z';q,uj) are 2D and 
3D interacting density response functions, respectively, 
(j)^^t(^z;q,uj) is the 2D Fourier transform of the exter- 
nal potential (/)'^^*(r, w), and v{z,z';q) is the 2D Fourier 
transform of the bare Coulomb interaction: 



«(z,z';g)=z;,e-'l^-^'l, 

with Vq — 271 /q. 

Combining Eqs. ^ and Q, we find 



(4) 



Sn2D{z;q,uj) ^ 5{z- Zd)Xeff{.q,^^)(l>{z]q,uj), (5) 
where 

X2D(g,w) 



Xeffiq,^) 



1 - X2D (g, [W{zd, Zd] q, uj) 



(6) 



W{zd, Zd'jqjLL') being the so-called screened interaction 
W{z, z';q,iL!) = v{z, z';q) + J dzi J dz2 



X v{z, zi;q) x?.d{zi,Z2] q, w) v{z2,z'; q), 



(7) 



and (j){z; q, uj) being the 2D Fourier transform of the total 
potential at z in the absence of the 2D sheet: 



iz;q,uj) 



dz" 



X X3D{z',z";q,Lu)]r'''{z";q,L^). 



5{z-z")+ / dz' v{z,z';q) 



(8) 



Eq. © suggests that the screening of the 3D subsystem 
can be incorporated into the description of the electron- 
density response at the 2D electron gas through the in- 
troduction of the effective density-response function of 
Eq. (jSJ, whose poles should correspond to 2D electron- 
density oscillations. 

Alternatively, we can focus on the 2D Fourier trans- 
form of the total potential at z in the presence of both 
2D and 3D subsystems: 



(/)(2:;g,a;) = (j)'''^*{z;q,uj) + J dz' v{z, z'; q) 

X [Sn2Diz';q,uj) + Sn3D{z';q,uj)], (9) 

which with the aid of Eqs. © and ((SJl can also be ex- 
pressed in the following way: 



(2) (j){z;q,uj) ^ (j){z;q,uj)+ J dz'W{z,z';q,uj)5n2D{z';q,uj). 

(10) 

Now we choose z = Zd, and using Eq. ij^l we write 



(j){zd;q,uj) = [l + W{zd,Zd;q,uj)xeff{q,i^)] 4>izd;q,u;), 

(11) 

which allows to introduce the effective inverse 2D dielec- 
tric function 



Jf{q,i^) = l + W{zd,Zd;q,uj)xeff{q,i^)- 



(12) 



Since our aim is to investigate the occurrence of long- 
wavelength (q —^ 0) collective excitations, we can rely on 
the random-phase approximation (RPA),^^ which is ex- 
act in the q limit. In this approximation, the 2D and 
3D interacting density-response functions are obtained as 
follows 



X2D(g,^) 



(13) 



and 



X3D{z,z';q,Lj) ^ xloiz^^'iQ^^) + / dz 



dzn 



X X3D(^7^i;'7'^)"(^l>^2;'7)X3D(22,z';q,tj), (14) 

where X2Dil'^) ^^'^ xtoi^'-^'^'i^^) represent their non- 
interacting counterparts. An explicit expression for the 
2D noninteracting density-response function X2Dil'^) 
was reported by Stern^ In order to derive explicit 
expressions for the 3D noninteracting density-response 
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function XsdC^? -2;'; q, w) one needs to rely on simple mod- 
els, such as the hydrodynamic or infinite-barrier model, 
but accurate numerical calculations have been carried 
ouii^^i^ from the knowledge of the cigenfunctions and 
eigenvalues of the Kohn-Sham hamiltonian of density- 
functional theory (DFT)ii& 

Combining Eqs. ®, 112), and (H^Jl, one finds the RPA 
effective 2D dielectric function 



(15) 



The longitudinal modes of the 2D subsystem, or plas- 
mons, are solutions of 



ee//(g,w) = 0. 



(16) 



In the absence of the 3D subsystem, the 3D screened 
interaction W{z, z';q,uj) reduces to the bare Coulomb 
interaction v{z,z';q), and the solution of Eq. (fT?^ leads 
at long wavelengths to the well-known square-root wave- 
vector dependence of the 2D plasmon energjii^ 



(17) 



qp and m being the 2D Fermi momentum and 2D effec- 
tive mass, respectively. The 2D Fermi velocity is simply 

Vp — qp/lTT'. 

In the presence of the 3D subsystem, the long- 
wavelength limit of the effective 2D dielectric function 
of Eq. H15|) is found to have two zeros. One zero cor- 
responds to a high-frequency {lh » vpq) oscillation in 
which 2D and 3D electrons oscillate in phase with one an- 
other. The other mode corresponds to a low-frequency 
acoustic oscillation in which both 2D and 3D electrons 
oscillate out of phase. 

At high frequencies, where lo >> vpq^ the long- 
wavelength limit of the 2D density-response function 
X2£)('?,w) is known to be 



1 Lui 



limX2D(9,w >> vpq) = 



(18) 



On the other hand, when the 2D sheet is located ei- 
ther far inside or far outside the metal surface, the 
long-wavelength limit of the 3D screened interaction 
W{zd, Zd] q,LLi) takes the form 



lim W{zd, Zd; q,LO » vpq) 

(j-»0 



(19) 



where cup^s represents either the bulk-plasmon frequency 
LOp = \/4:Trn or the conventional surface-plasmon en- 
ergy ujs — tj-'p/V^,^'^ depending on whether the 2D sheet 
is located inside or outside the solid. Introduction of 
Eqs. fll^ and f]ll|l into Eqs. l(T5|l and l(TH|l yields a high- 
frequency mode at 



: CO. 



+ ^lo- 



rn 



At low frequencies, we seek for an acoustic 2D plasmon 
energy that in the long-wavelength limit takes the form 



a vp q. 



(21) 



A careful analysis of the 2D density-response function 
X2£i(9, w) and the 3D screened interaction W{zd, Zd; q, oj) 
shows that at to = avp q the long-wavelength limits of 
these quantities take the form 



limX2D(5',ai^i=^9) = - 



- 1 



and 



lim W{zd, Zd; q, avpq) = I{zd) 



(22) 



(23) 



An inspection of Eqs. (|15|l . (|22|l . and (|23|l indicates that 
for a low-energy acoustic oscillation to occur the quantity 
I{zd) must be different from zero. In that case, the long- 
wavelength limit of the effective 2D dielectric function 
of Eq. (|15|l has indeed a zero corresponding to a low- 
frequency oscillation of energy given by Eq. (|21|l with 



a = \ 1 



7r[7r + 2/(zrf)]' 



(24) 



In the following, we investigate the impact of the 3D 
screening on the actual wave-vector dependence of the 
low-energy 2D collective excitation. We first consider the 
two limiting cases in which the 2D sheet is located far in- 
side and far outside the metal surface, and we then carry 
out self-consistent calculations of the 3D screened inter- 
action W{z, z';q,uj), which will allow us to obtain plas- 
mon dispersions for arbitrary locations of the 2D sheet. 



A. 2D sheet far inside the metal surface 

In the case of a 2D sheet that is located far inside the 
metal surface, the 3D subsystem can safely be assumed 
to exhibit translational invariance in all directions, i.e., 
the screened interaction W{zd, Zd] q, uj) entering Eq. H15() 
can be easily obtained from the knowledge of the inter- 
acting density-response function X3D{k,u!) of a uniform 
3D electron g follows 



W{zd,zd;q,uj) 



dfh -1 



(fc,w), 



(25) 



where k = \/ q^ + q^ is the magnitude of a 3D wave vec- 
tor and e^^{k,uj) is the inverse dielectric function of a 
uniform 3D electron gas: 



e^^ik, uj)^l + — X3D{k, w) 



In the RPA, 



esoik, w) = 1 - xloik, t^), 



(26) 



(27) 



x'^pi{k,uj) being the noninteracting density-response 
function first obtained by Lindhardii^ 



1. Local 3D response 

If one characterizes the 3D uniform electron gas by a 
local dielectric function 63/3 (w), then Eq. (|25|l yields 

W'"-'^\za,Zd-q,Lo)^v^e-^{Lo). (28) 

In a 3D gas of free electrons, e3_D(w) takes the Drude form 



e3D(w) = 1 2' 



which yields 



(29) 



(30) 



This means that in a local picture of the 3D response 
the characteristic collective oscillations of the 2D electron 
gas would be completely screened by the sorrounding 3D 
substrate and no low-energy acoustic mode would exist ^ 



2. Hydrodynamic 3D response 

Dispersion effects of the 3D subsystem can be incorpo- 
rated approximately in a hydrodynamic model. In this 
approximation, the dielectric function of a 3D uniform 
electron gas is found to bei^ 



/32 p 



(31) 



where P = ^JTJ?,kF represents the speed of propagation 
of hydrodynamic disturbances in the electron system, 
and fci? is the 3D Fermi momentum. 

Introducing Eq. (|31|l into Eq. H25() . one finds 

lim W{zd, Zd] g, avpq) = 27r/3/wp, (32) 

which yields the following simple expression for the 
acoustic coefhcient of Eq. (OH : 



a=^ 1 



I + API Up' 



(33) 
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FIG. 1: Effective dielectric function of a 2D sheet that is 
located far inside the metal surface, as obtained from Eq. I15II 
with (a) q = 0.01 and (b) q = 0.1. The real and the imaginary 
parts of e^fj{q,u)) are represented by thick and thin solid 
lines, respectively. The dotted line represents the effective 2D 
energy- loss function Im ejj?^(g, tj)] . The vertical dashed 
line represents the upper edge uju ~ vpq + q^ /{2m) of the 
2D e-h pair continuum, where 2D e-h pairs can be excited. 
2D and 3D electron densities have been taken to be those 
corresponding to the Wigner radii r1^ = 3.14 and = 1.87, 
respectively. The 2D effective mass has been taken to be 
m = 1. 



3. Full 3D response 

We have carried out numerical calculations of the RPA 
effective dielectric function of Eq. H15|) , by using the full 
X2d(9i'^) ^ii'i XsoikT^^) density-response functions, and 
choosing the electron-density parameters r^^ = 3.14 and 
rf-^ — 1.87 corresponding to the (0001) surface of Be.^^ 

The results we have obtained with q = 0.01a(^^ and 
q = O.log ^ are displayed in Figs. ^ and respec- 
tively. We observe that at energies below the upper edge 
= vpq + q^ /{"im) (vertical dashed line) of the 2D e- 
h pair continuum (where 2D e-h pairs can be excited) 



the real part of the effective dielectric function is nearly 
constant and the imaginary part is large, as would occur 
in the absence of the 3D susbtrate. At energies above 
cij„, momentum and energy conservation prevents 2D e- 
h pairs from being produced, and Imee//((7, w) is very 
small. 

Collective excitations are related to a zero of 
Reee//((?, w) in a region where Imeef f{q,u:) is small and 
lead, therefore, to a maximum in the energy-loss func- 
tion lTa[—e~^^{q,u})]M' In the absence of the 3D sub- 
strate, a 2D plasmon would occur at ix}2D = 1.22 eV for 
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FIG. 2: Stars represent the q coefBcient of the acoustic- 
plasmon energy uj = avpq versus the 3D Wigner radius, as 
obtained from Eq. 11611 in the long- wavelength limit. These 
results are found to be insensitive to the 2D Wigner radius 
r^^ . The solid line represents the prediction of Eq. 12411 . as 
obtained with the full RPA value of I{zd — > —<x). The dotted 
line represents the hydrodynamic prediction of Eq. 1331 . 



q — O.OIoq ^ and uj2d = 3.99 eV for q — O.Ioq ^. However, 
Fig. n shows that in the presence of the 3D substrate 
a well-defined low-energy acoustic plasmon occurs, the 
sound velocity being just over the 2D Fermi velocity vp- 
The small width of the plasmon peak is entirely due to 
plasmon decay into e-h pairs of the 3D substrate. 

We have carried out calculations of the effective 2D di- 
electric function of Eq. (|15|l for a variety of 2D and 3D 
electron densities, and we have found that a well-defined 
acoustic plasmon of energy co = avp q always present 
for 2D wave vectors up to a maximum value of g ^ 
where the acoustic-plasmon dispersion merges with lUu- 
The coefficient a that we have obtained from the zeros 
in Eq. (|16|l is represented by stars in Fig. [3 versus the 
3D Wigner radius r^^, together with the prediction of 
Eq. H24|l as obtained with the computed RPA value of 
1(^(1 — * — oo) (solid line) and the hydrodynamic predic- 
tion of Eq. (1231 (dotted line). Fig. |3 shows that Eq. ^ 
is a good representation of the linear dispersion of this 
low-energy plasmon, especially at the highest 3D electron 
densities. Fig. |2 also shows that for low electron densi- 
ties the hydrodynamic prediction is too small, which is 
due to the fact that at low densities the long-wavelength 
limit of the 3D screened interaction is underestimated in 
this approximation. 



B. 2D sheet far outside the metal surface 

In the case of a 2D sheet that is located far outside the 
metal surface, where the 3D electron density is negligible, 
the 3D screened interaction of Eq. ^ at z = z' = Zd 



takes the form 

W{zd,zd;q,uj) = vg[l- e~2 1 g{q^ ^)] , (34) 

where g{q, lo) is the so-called surface-response function of 
the 3D subsystem 



dzi / dz2e«('^i+^^'x3D(2i,22;'?,w)- 



(35) 



Local 3D response 



In the simplest possible model of a metal surface, one 
characterizes the 3D substrate at z < by a local dielec- 
tric function which jumps discontinuously at the surface 
from e^oi^) inside the metal (z < 0) to zero outside 
{z > 0). Witin this modcl,^^ 



e3D(^) ~ 1 



(36) 



which is precisely the long- wavelength limit of the actual 
surface-response function. 

At low frequencies, where e^ni^^) is large [see Eq. H29|l ] 
and g^^^^'iq^uj) ^ 1, Eq. (^3 yields 



lim W^""'"'' 



{zd,Zd; q^avpq) ^ 4:Tr Zd. (37) 
Introducing Eq. H37|) into Eq. H24|l . one finds 



a 




(38) 



For large values of the distance Zd between the 2D sheet 
and the metal surface, one can write 

a « V2l^, (39) 

which is the result first obtained by Chaplik^ by using 
the Drude-like 2D density-response function of Eq. 



2. Nonlocal 3D response 

An inspection of Eq. (|34|l shows that the 
long-wavelength limit of the screened interaction 
W{zd, Zd;q,io) is dictated not only by the local {q = 0) 
surface-response function g'°'^°'(g,w) but also by the 
leading correction in q of the actual nonlocal g(q,Lu). 
Feibelman showed that up to first order in an expansion 
in powers of q, the surface-response function of a jellium 
surface can be written asS^ 



5(9, ^) 



1 + 2qd±{iL;] 



e3D(w) 



£315 (w) + 1 



0{q% 
(40) 



which at low frequencies yields 

5(g,a;)«l + 2gdi(0). 



(41) 



The frequency-dependent d±{uj) function occurring in 
Eq. H40|l represents the centroid of the induced 3D charge 
density, which in the static limit (cu = 0) reduces to the 
image plane of an external point charge. 

Using Eq. ||1T|| . we find the actual long- wavelength 
limit of Eq. 



lim W{zd, Zd] q, avpq) = 4 tt [z^ - 
which combined with Eq. H24() yields 



d±(0)] 



a 




(42) 



(43) 



This shows that the acoustic-plasmon sound velocity de- 
rived from the local model [see Eq. H38() ] remains un- 
changed, as long as is replaced by the coordinate of 
the 2D sheet relative to the position of the image plane. 



3. Full 3D response 

In order to compute the full RPA surface-response 
function of Eq. H35|l . we follow the method described 
in Ref. Q for a jellium slab. We first assume that 
the 3D electron density vanishes at a distance zq from 
either jellium edge,^^ and compute the noninteracting 
density-response function Xsni^^ -^'i 9i ^) from the knowl- 
edge of the self-consistent Kohn-Sham wavefunctions and 
energies of DFTji& which we obtain in the local-density 
approximation (LDA)i2& We then introduce a double- 
cosine Fourier representation for both the noninteract- 
ing and the interacting density-response functions, and 
find explicit expressions for the surface-response func- 
tion in terms of the Fourier coefhcients of the density- 
response function.— To ensure that our slab calcula- 
tions are a faithful representation of the actual surface- 
response function of a semiinfinite 3D system, we follow 
the extrapolation procedure described in Ref. '28'. 

We have carried out numerical calculations of the ef- 
fective dielectric function of Eq. H15|l . by using the full 
2D noninteracting density-response function, X2Dil^^)^ 
and the self-consistent RPA surface-response function, 
g{q,u!), with electron-density parameters r^^ — 3.14 and 
r^D = 1.87 corresponding to Be(OOOl). 

The results we have obtained for a 2D sheet located at 
Zd = Af are displayed in Figs.l^L (with q = O.Ola,^^) and 
^jp (with q = O.Igq^), \f = 2'K/kp being the 3D Fermi 
wavelength. Fig. O clearly shows that in the presence 
of the 3D substrate a well-defined low-energy acoustic 
plasmon occurs, the sound velocity being close to that 
predicted by Eq. (gSJ with d±{0) = 0.2 Af (vertical long- 
dashed lines). The actual plasmon energy is smaller than 
predicted by Eq. H43|l . especially at the shortest wave- 
lengths (q = O.Icq ^), simply due to the bending of the 
plasmon dispersion as a function of q (see Fig. below). 



3 





1 1 1 1 


(a) 


: J 


















0.1 0.2 0.3 0.4 

CO(eV) 



3 



1 


1 1 1 


(b) 




1 


1 . 1 , 1 



1 2 3 4 5 

ro(eV) 

FIG. 3: As in Fig.Q but now for a 2D sheet that is located at 
one 3D Fermi wavelength outside the metal surface {zd — Xf)- 
The long-dashed vertical lines here represent the plasmon- 
energy uj — avp q predicted by Eq. ^ with d± = 0.2 Af. 
For real frequencies, a 2D sheet that is located a.t Zd ~ Xf 
exhibits a plasmon peak that at q — 0.01 is extremely 
sharp (as — > oo the plasmon peak becomes a delta func- 
tion); hence, in the calculations presented in this figure we 
have replaced the energy tJ by a complex quantity u + irj with 
(a) r; = 0.05 eV for q = 0.01 a^^ and (b) 77 = for q = 0.1 a^^. 



C. 2D sheet at an arbitrary location 



Hydrodynamic 3D response 



An explicit expression for the screened interaction 
W{z,z';q,tj) of Eq. (Q can be obtained in a hydrody- 
namic model in which the 3D electron density is assumed 
to change abruptly at the surface from n inside the metal 
to zero outside. After writing and linearizing the ba- 
sic hydrodynamic equations, i.e., the continuity and the 



Bernouilli equation, we find 

r (1 _ e2-<^-p//5) 2^ < 
limW{zd,Zd;q,avpq) = < 

q — s-0 I 

L 4 TT Zd> 

(44) 

which combined with Eq. H24() yields an explicit expres- 
sion for the acoustic coefficient a. We note that in a local 
description of the electronic response of the solid surface 
(/3 = 0) the 3D screened interaction W{zd, Zd] q, avpq) is 
zero inside the solid {zd < 0) and An Zd outside {zd > 0). 
This shows that in the 2D long-wavelength limit {q 0) 
the nonlocality of the 3D response is only present inside 
the solid {zd < 0), where finite values of the 3D momen- 
tum k are possible. 

Alternatively, the screened interaction W{z, z'; q,uj) 
can be obtained within a specular-reflection model 
(SRM)22 Qj-^ equivalently, a classical infinite-barrier 
model (CIBM)^°'^^ of the surface, which have the virtue 
of describing the 3D screened interaction in terms of the 
bulk dielectric function e^nik, uj) of a 3D uniform (and in- 
finite) electron gas (see Appendix 0. If this bulk dielec- 
tric function is chosen to be the hydrodynamic dielectric 
function of Eq. (jSIJ , then these models yield Eq. . A 
more accurate description of the bulk dielectric function 
e3D{k,uj) yields a result that still coincides with that of 
Eq. 144|) outside the surface {zd > 0), though small dif- 
ferences may arise at z^j < 0. 

When the 2D sheet is located far inside the metal 
{zd << 0), Eq. (|m yields the hydrodynamic asymptotic 
behaviour dictated by Eq. (|32|l . and the SRM combined 
with the RPA bulk dielectric function yields the correct 
RPA asymptotic behaviour. However, these hydrody- 
namic and specular-reflection models, which are both 
based on the assumption that the 3D electron density 
drops abruptly to zero at the surface, fail to reproduce 
the correct asymptotic behaviour outside the surface [see 
Eq. H42I) ]. This is due to the fact that the leading cor- 
rection in q of the surface-response function g{q,uj) is 
governed by the spill out of the electron density into the 
vacuum, which is not present in these models. 



2. Full 3D response 

For an arbitrary location of the 2D sheet we need to 
compute the full screened interaction W{zd, Zd^qjUj) of 
Eq. ((TJ. To calculate this quantity we consider a jellium 
slab, as we did to obtain the surface-response function 
g(q,uj), and we find explicit expressions in terms of the 
Fourier coefficients of the interacting density-response 
function, which we compute in the RPA [see Eq. (|14|l ] 
from the knowledge of the LDA eigenvalues and eigen- 
functions of the Kohn-Sham hamiltonian of DFT. 

In Fig. ^ the long-wavelength limit I{zd) of the 
screened interaction W{zd, Zd;q,avFq) [see Eq. 1^5)1 ] is 
displayed versus Zd, as obtained with rf^ = 1.87 from 
our full self-consistent RPA calculations (thick solid line) 
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FIG. 4: Long-wavelength limit I{zd) of the screened interac- 
tion W{zd, Zd', q, avpq)- The thick solid line represents the full 
self-consistent RPA calculation. The results obtained from 
Eq. (EU with dx(0) = 0.2 Af and from Eq. ^ are repre- 
sented by dotted and thin solid lines, respectively. The hor- 
izontal dashed line represents the result obtained from the 
RPA bulk screened interaction of Eq. II25II . The 3D Wigner 
radius has been taken to be rf^ — 1.87. 



and from the hydrodynamic Eq. H44I) (thin solid line). 
Far inside the solid, our full calculation is close to the 
hydrodynamic prediction (see also Fig. [SJ and coincides 
with the result one obtains from the bulk screened in- 
teraction of Eq. (|25|) (horizontal dashed line). Near the 
surface, our full calculation considerably deviates from 
the hydrodynamic prediction and converges far outside 
the solid with the asymptotic curve dictated by Eq. H42|l 
with d_L(0) = 0.2 Af (dotted line).^ 

At this point, it is interesting to note that within a lo- 
cal picture of the 3D response the long-wavelength I{zd) 
screened interaction would be zero for all locations of 
the 2D sheet inside the metal {zd < 0), showing that 
the screening of 2D electron-density oscillations would 
be complete and no acoustic surface plasmon would oc- 
cur. It is precisely the nonlocality of the 3D response 
(finite values of the 3D momentum k are still present in 
the 2D long-wavelength limit) which provides incomplete 
screening and allows, therefore, the formation of acoustic 
surface plasmons in the interior of the solid. We also note 
that within a simple nonlocal picture of the 3D response, 
such as the hydrodynamic and specular-reflection mod- 
els described above, the screening of 2D electron-density 
oscillations would still be complete at the jellium edge 
{zd — 0). Hence, in the real situation where the 2D 
surface-state band is located very near the jellium edge 
the occurence of acoustic surface plasmons is originated 
by a combination of the nonlocality of the 3D response 
and the spill out of the 3D electron density into the vac- 
uum. 

Figs. and \^ exhibit the results we have obtained 
for the effective dielectric function of Eq. (|15|l [with 
q = 0.01 a(7^ (Fig.ISi) and q = 0.1 a^^ (Fig. Eb)] by us- 
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FIG. 5: As in Fig. Q but now for a 2D sheet that is lo- 
cated at the jellium edge {zd — 0). Also shown is the effec- 
tive 2D energy-loss function Im for = — Af, 
Zd = and Zd = Xf (dotted lines). The open circles rep- 
resent the effective 2D energy-loss function Im ejj?^(q, a;)] 
obtained from the limiting Eq. 1251 appropriate for a 2D 
sheet far inside the metal and from the limiting Eq. 1341 with 
Zd = Xf appropriate for a 2D sheet far outside the metal. 
These calculations are found to coincide with the full calcula- 
tions for Zd — —Xf and Zd = Af, respectively. As in Fig.l^Ja), 
the calculations presented here for Zd = Xf and q = 0.01 a^^ 
have been carried out by replacing the energy by a complex 
quantity u; + irj with rj — 0.05 eV. All remaining calculations 
have been carried out for real frequencies, i.e., with 77 = 0. 



ing the full 2D noninteracting density-response function, 
X^^jC^jw), and the self-consistent RPA screened inter- 
action, W{zd, Zd', q,Lij), with electron-density parameters 
r^^ = 3.14 and rf^ = 1.87 corresponding to Be(OOOl). 
In these figures the 2D sheet has been taken to be lo- 
cated at = 0, as approximately occurs with the quasi- 
2D surface-state band in Be(OOOl). For comparison, also 
shown in these figures are the results we have obtained 
for the energy-loss function when the 2D sheet is located 



FIG. 6: The open circles represent the sound velocity Vs 
{uj — Vs q) of the low-energy acoustic plasmon that is visible 
in Fig. |K| versus the location Zd of the 2D sheet with respect 
to the jellium edge. The horizontal short-dashed line repre- 
sents the result we have obtained from the limiting Eq. 12511 
appropriate for a 2D sheet far inside the metal. The dotted 
line represents the result we have obtained from the limiting 
Eq. 14311 with d± — 0.2 Xf, which is appropriate for a 2D 
sheet far outside the metal. The long-wavelength limit vf of 
the upper edge uju/q of the 2D e-h pair continuum is repre- 
sented by an horizonal long-dashed line. The thick and thin 
solid lines represent the results obtained from Eq. 124|l by us- 
ing the actual RPA I{zd) and the hydrodynamic Eq. 144L 
respectively. 2D and 3D electron densities have been taken 
to be those corresponding to the Wigner radii = 3.14 
and rl^ — 1.87, respectively. The 2D effective mass has been 
taken to be m = 1. 



inside the metal at Zd — — Af and outside the metal at 
Zd — Af/2 and Zd = Xp. 

An inspection of Fig. shows that (i) the results we 
have obtained for Zd — — Af and Zd = Xp are exactly 
reproduced by the limiting Eqs. H25|l and (|34|l appropriate 
for a 2D sheet far inside and far outside the metal surface, 
respectively; and (ii) in the actual situation where the 
2D surface-state band is located very near the jellium 
positive background edge {zd = 0), a well-defined low- 
energy acoustic plasmon occurs, the sound velocity being 
very close to the limiting case of a 2D sheet far inside 
the metal surface and being, therefore, just above cuu- 
This is in agreement with the recent prediction that in 
a real metal surface where a partially occupied quasi-2D 
surface-state band coexists in the same region of space 
with the underlying 3D continuum an acoustic surface 
plasmon should occur at energies just above the upper 
edge of the 2D e-h pair continuum.^ 

The sound velocity Vs (w = Vs q) of the acoustic plas- 
mon that is visible in Fig. is displayed in Fig. versus 
the location Zd of the 2D sheet relative to the jellium 
edge, as obtained from our full RPA self-consistent cal- 
culation of the effective 2D dielectric function of Eq. H15|) 
(open circles), together with the sound velocity Vg = avp 
obtained from Eq. (021 with d^(0) = 0.2 Af (dotted line). 
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FIG. 7: Dispersion of the acoustic plasmon occurring in a 2D 
sheet that is taken to be located far inside the solid (thick 
dotted line), ai Zd = (open circles), at Zd ~ Xp (solid line), 
and infinitely far outside the metal (solid circles). The thick 
dashed line represents the upper edge uju = vpq + /(2m) of 
the 2D e-h pair continuum. The thin dotted line represents 
the 2D plasmon energy lo2d dictated by Eq. 1)17^ . which is 
accurate at long wavelengths (g 0). The thin dashed line 
represents the 2D plasmon energy yj u)2d + 3'^F'i'^/4 that is 
obtained after an expansion of 'x^oi'i^'^) iii powers of VFqjuj. 
2D and 3D electron densities have been taken to be those 
corresponding to the Wigner radii r^^ = 3.14 and = 1.87, 
respectively. The 2D effective mass has been taken to be 
m = 1. 



When the 2D sheet is located inside the metal surface, 
the sound velocity nicely converges with the RPA bulk 
calculation from Eq. H25I) (horizontal short-dashed line). 
When the 2D sheet is located outside the metal surface, 
the sound velocity converges with the limiting value a vp 
obtained from Eq. and dj_{0) = 0.2 Xp- For com- 
parison, also shown in this figure is the result we have 
obtained from Eq. by using the actual RPA I{zd) 
screened interaction (thick solid line) and from the hy- 
drodynamic Eq. 1)44(1 . These calculations clearly show 
that Eq. H24I) accurately reproduces the dispersion of 
acoustic surface plasmons, as long as the long- wavelength 
limit I{zd) of the screened interaction is described self- 
consistently with full inclusion of the electronic selvage 
structure at the surface. 

The sound velocity of Fig. (open circles) has been 
obtained from the effective 2D dielectric function at very 
low 2D momenta, where the energy of the acoustic plas- 
mon is linear in q. The behaviour of this plasmon en- 
ergy as a function of the 2D momentum q is displayed in 
Fig. [7| with the 2D sheet chosen to be located far inside 
the solid (thick dotted line), at = (open circles), at 
Zd = (solid line), and infinitely far outside the solid 
(solid circles) . The upper edge of the 2D e-h pair contin- 



uum is represented by a thick dashed line, showing that in 
the real situation where the 2D sheet is located near the 
jellium edge the energy of the acoustic surface plasmon 
(open circles) is just outside the 2D e-h pair continuum 
for all momenta under study. 

III. SUMMARY AND CONCLUSIONS 

The partially occupied band of Shockley surface states 
in a variety of metal surfaces is known to form a quasi-2D 
electron gas that is immersed in a semiinfinite 3D gas of 
valence electrons. In order to describe the impact of the 
dynamical screening of the semi-infinite 3D continuum 
on the electronic excitations at the 2D electron gas of 
Shockley surface states, we have presented a model in 
which the dynamical screening of 3D valence electrons is 
incorporated through the introduction of an effective 2D 
dielectric function. 

We have considered the two limiting cases in which 
the 2D sheet is located far inside and far outside the 
metal surface. In both cases, the dynamical screening 
of the valence electrons in the metal is found to change 
the 2D plasmon energy from its characteristic square- 
root behaviour to a linear dispersion, the sound velocity 
being proportional to the Fermi momentum of the 2D 
gas. As this collective oscillation occurs in a region of 2D 
momentum space where 2D e-h pairs cannot be produced, 
this is a well-defined acoustic plasmon. The finite width 
of the plasmon peak is due to a small probability for the 
plasmon to decay into e-h pairs of the 3D substrate. 

We have shown explicitly that when the 2D sheet co- 
exists in the same region of space with the underlying 
3D continuum the origin of acoustic surface plasmons, 
which have been overlooked over the years, is dictated by 
a combination of the nonlocality of the 3D response and 
the spill out of the 3D electron density into the vacuum, 
both providing incomplete screening of the 2D electron- 
density oscillations. 

We have carried out self-consistent DFT calculations 
of the dynamical density-response function of the 3D sys- 
tem of valence electrons, and we have found that a well- 
defined acoustic plasmon exists for all possible locations 
of the 2D sheet relative to the metal surface. The en- 
ergy dispersion of this acoustic surface plasmon is slightly 
higher than the energy of the collective excitation that 
has recently been predicted to exist at real metal sur- 
faces where a quasi- 2D surface-state band coexists with 
the underlying 3D continuum^ Small differences between 
the plasmon energies obtained here and those reported 
previously*^" are due to the absence in the present model 
of transitions between 2D and 3D statesi^ 
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called classical infinite-barrier model (CIBM) of a jellium 
surface rSiSi one finds 



APPENDIX A: SPECULAR-REFLECTION 
MODEL OF THE 3D RESPONSE 



Either by assuming that electrons are specularly re- 
flected at the surface (SRM)^^ or by invoking the so- 



W{zd, Zd;q,uj) = Vq 



l_e-2<?-. [1 -e,(0;^)]/ [1 + 6,(0; c^)], 



Zd>0 



es{0; q,Lo) + es{2zd;q,uj) - 2el[zd]q,uj)/ [es(0;g,w) -I- 1] , elsewhere 

I 



(Al) 



where k = v?+^ ^ momentum, and 

dq 



es{z;q,uj) 



q 



^e^'^^e3^(fc,c.), 



(A2) 



^3D{q,'-^) being the dielectric function of a uniform (and 
infinite) 3D electron gas. 

If the 3D dielectric function e3£)((7, oj) is chosen to be 
the hydrodynamic dielectric function of Eq. H31() . then 



one finds 



es{z]q,uj) = 



9 9 
tU^ — LU^ 



P^lq 



-q\zd\ 



which in combination with Eq. IjAip yields Eq. (|44|l . 



(A3) 
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